Considerations about the Aharonov- Anandan Phase for Time Independent 

Hamiltonians 



o 
o 

(N 



Oh' 
I 

> ■ 

■ 

a ■ 

=3 : 



O 
m 
o 



oo 

o 



P.-L. Giscard 

B2 Institute, Department of Physics and College of Optical Sciences 
The University of Arizona, Tucson, Arizona 85721 
(Dated: January 3, 2009) 

We present a method for calculating the Aharonov-Anandan phase for time-independent Hamil- 
tonians that avoids the calculation of evolution operators. We compare the generic method used to 
calculate the Aharonov-Anandan phase with the method proposed here through four examples; a 
spin-i particle in a constant magnetic field, an arbitrary infinite-sized Hamiltonian with two known 
eigenvalues, a Fabry-Perot cavity with one movable mirror and a three mirrors cavity with a slightly 
transmissive movable middle mirror. 
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I. INTRODUCTION 

The geometric phase, which was first recognized to 
arise from quantum systems undergoing a cyclic mo- 
tion by M. V. Berry as been studied in much detail, 
both from the mathematical point of view and for vari- 
ous physical systems. The topological nature of Berry's 
phase and the Aharonov-Anandan phase has been inves- 
tigated since its discovery, yielding deep insights in the 
cyclic motion of quantum systems 0, M, 0|- De- 
spite a good understanding of these geometric phases Q , 
[7| and their use in various branches of physics, from 
quantum phase transitions Q to quantum computation 
(see for example 0], few studies have focused on 

algebraic methods to calculate them [llj, [l^]. In fact 
the generic method still in use to calculate geometric 
phases was given in the pioneering article by Aharonov 
and Anandan [l3j . 

The main goal of this paper is to propose an alterna- 
tive method for calculating geometric phases arising from 
time independent Hamiltonians that avoids the explicit 
calculation of the evolution operator. We also derive re- 
sults concerning the condition for an Hamiltonian to give 
rise to a cyclic motion of a state and on the possible re- 
sulting total phases (see Sec lH B[) . Sec lHII then gives four 
examples of calculation of the geometric phase, showing 
that the results coincide with those obtained using the 
generic way to determine it. This also enables a compar- 
ison of the complexity involved in each method. Finally, 
in the last example we calculate the period and the to- 
tal phase of a three-mirror cavity, using both the generic 
method and the new method proposed here. 



II. CALCULATING THE GEOMETRIC PHASE 
A. A simplified expression for the geometric phase 

Let us first review the technique used to calculated the 
geometric phase arising from a cyclic motion of a quan- 
tum system. A basic idea for the calculation of the geo- 



metric phase is to calculate the dynamical phase (ipDyn) 
and to subtract it from the total phase accumulated dur- 
ing one cyclic evolution of the system under considera- 
tion. In the following we only study the calculation of 
the Aharonov-Anandan phase (AA-phase). 

Let T be the period of this evolution and the phase 
difference between the initial state |^'(0)) and |^'(r)), 
i.e. |*(t)) = e*'^|^'(0)). Now let f{t) be a continuous 

and let |*(i)) = 



function of time so that /(r) — /(O) 



that the geometric phase 7 is given by 



\E'(t)). Aharonov and Anandan have proven in [ij 
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(1) 



For practical reasons, one uses instead the equivalent fol- 
lowing formula that involves explicitly 



J^(b + l {^{t)\^it))dt = ^-ipDy„, (2) 

Jo 

where one has to calculate \'^{t)) and then |4'(t)) from 
the Hamiltonian, using the evolution operator U{t). If for 
any time t and t' , [H{t), H{t')] ^ 0, then U{t) = p-^Ht/n^ 
However for complicated systems, calculating U (t) is not 
trivial [IBl and the methods used for disentangling 
it are generally cumbersome [iGj , [l3| • Furthermore once 
a disentangled form of U has been found, its action on a 
given initial state |^'(0)) may lead to involved expressions 
for |4'(t)). But as we show later on, these explicit calcu- 
lations can be avoided if the Hamiltonian can be written 
as H'{t) — H{t) + \a{t) where a{t) is a function of time 
(not an operator) and H{t) obeys [H{t),H{t')] = 0. 

As the geometric phase arising from a given Hamilto- 
nian is gauge invariant [l3]), the geometric phases 
arising from H'{t) and H{t) are the same. Using 
i|*(t)) = h-^H{t)\^{t)) and \^(t)) = J7(i)|*(0)) the 
geometric phase is therefore 

-f = cb + h-' f {'f{Q)\U^{t)H{t)U{t)\'i>{0))dt. (3) 



For Hamiltonians that commute with themselves at 
any time, [U{t'),H{t)] = [W{t'),H{t)] = and since 
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W{t)U{t) = U{t)W{t) = 1, the geometric phase reduces 
to 

7 = + ?i"'/ (*(0)|i/(OI*(0))di. (4) 
Jo 

This expression can be further simphfied in the case of 
time independent Hamihonians to 

7 = 0+^(vI/(O)|/J|*(O)). (5) 

Since the system evolves according to the Schrodinger 
equation, that is to say at constant energy, and is driven 
by a time independent Hamiltonian, r and are the same 
for any state located on the closed loop formed by the 
cyclic time-evolution in the projective Hilbert space. In 
other terms, starting the cyclic evolution from any point 
of the loop always yields the same geometric phase. 

B. Finding the Period and Total Phase of the 
Cyclic Motion 

1. The Period r 

In the following we only consider time-independent 
Hamiltonians. The period r and the phase are gen- 
erally derived from physical considerations and a careful 
study of the evolution operator U{t). One can however 
avoid such a study by working only with the Hamilto- 
nian. Let B = {|0fc)} be a basis in which H is diag- 
onal and A = {Afc} the corresponding set of its eigen- 
values. Let i?* C i? be the smallest set of eigenvectors 
needed to decompose a state \^) on B and let A$ C A 
be the corresponding set of eigenvalues. Finally, let 
AE'* be the set of non-zero energy spacings in A^,, i.e. 

Then, we show in the Appendix [XI that the period r of 
the cyclic motion for 1 4') is proportional to the least com- 
mon multiple (LCM) of the inverse of AE^ . Rigorously 
we find 

T = 2TTh LCM {AE^^) . (6) 

A simple example of application of Eq.([5]) is the cal- 
culation of the period of evolution of a coherent state 
driven by a free field Hamiltonian (e.g. H = hioa^a). 
Eq.® gives 

T = 2Tih LCM {{hwn)-^) ^ 2tt/uj, (7) 

where we have used the fact that the LCM of all 1/n 
(with n a non-zero integer) is 1. 



2. Conditions of C'yclicality 

The form of the period r of Eq.® provides a means 
to test if a given Hamiltonian will yield a cyclic evolution 



of a given initial state. Indeed, if the system initially in 
state 1 4') effectively undergoes a cyclic motion, then r 
must be finite, i.e. the least common multiple involved 
in Eq.® must be finite. 

Therefore, there must be no two different elements of 
AE^^ that are incommensurable. In the case where Aij, 
contains more than two different eigenvalues [Tsl; this 
means that if one of the eigenvalues in A^ is irrational 
(say equal to k), then the only way H could yield a cyclic 
evolution of the state j^") is that all the eigenvalues in A>i, 
are rational multiples of k. 

Another condition for infinite Hamiltonians is that the 
set AE^^ must be bounded in E. This means that there 
must not be infinitely close eigenvalues in A^,. 

Remarkably, these conditions are always fulfilled if A^j, 
contains exactly two different eigenvalues. 

3. Total Phase 4> 

In the same way than for r, we show in Appendix |^ 
that the total phase accumulated by an initial state 
after one cycle of evolution is given by 

= 27r [n - A LCM {AE^^)] , (8) 

where n is an integer that depends on A G A.^. This 
expression being valid for any eigenvalue of A of A^, as 
soon as € A^, then the total phase must be = 2tt. 

In the case where A$ contains more than two different 
eigenvalues, Eq.® also yields a more general condition 
on <f>. Indeed, as we have seen before, the least common 
multiple involved in the above equation must be finite 
for the system to exhibit a cyclic motion of the state 
I*). Thus in Eq.®, A LCM (A^;^^) is always rational 
as soon as the system exhibits a cyclic motion, so that 

(f> 

— is rational. (9) 

TT 

This condition does not holds if A^, contains exactly 
two different eigenvalues (in which case (^/tt can be ir- 
rational). 

4- A Method for Time Independent Hamiltonians 

Ideally the geometric phase is just given by Eq.®. 
Nevertheless, as we see from Eq.® and Eq.®, r and 
<j) require a complete knowledge of the spectrum A.p in- 
volved in the decomposition of j^E') on the base B. It is 
however generally possible to use or impose relations on 
H so that the least common factor is easier to calculate 
(such assum ptio ns are often made in the generic method, 
see e.g. Ref.[l3). 

In the case where not all the eigenvalues are known, 
and there is no relation that one can or wants to im- 
pose on H to calculate the AA-phase, it is still possible 
to partially know it from only one non-zero element of 
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A^i. If the only known eigenvalue is zero, we have no 
information on r but = 27r is completely determined. 

Let A G Aij,. As shown in Appendix lAl there exists an 
integer n such that 



2mT. 



(10) 



Performing the gauge transformation H' ~ H + 
which leaves the AA-phase unchanged, Eq. (fTO| becomes 
h~^\T = —(j> and Eq.(I5]) transforms to 

7 = </-[l-A"i(*(0)|i7|*(0))] +27m, (11) 

which is consistent with the gauge invariance of the geo- 
metric phase. Note that the total phase </) in that equa- 
tion is unknown. It is also possible to recast the expres- 
sion of the AA-phase in term of an unknown period r 

as 



7 



[(*(0)|iJ|^'(0)) -A] + 27rn. 



(12) 



These relations hold for any element of A*. Only the 
total phase or the period remains unknown. However, if 
another element of Aip is known, it is possible to constrain 
the possible total phases or periods using Eq.((8]) for the 
two eigenvalues, see Sec lIIIBl In the opposite way, if the 
total phase or the period is known from physical consid- 
erations or if G Aij, , the geometric phase is given by the 
knowledge of only one non-zero eigenvalue and it is possi- 
ble to constrain the unknown part of the spectrum of the 
considered Hamiltonian. These procedures are detailed 
in Seclml 



III. SOME EXAMPLES 

In this section, we illustrate the calculation of the 
Aharonov-Anandan phase for various physical systems. 
We assume that the considered systems effectively un- 
dergo cyclic motions when starting in the proposed states 
|*(0)). 



A. Spin-i 

Consider first the precession of a spin-i around a con- 
stant magnetic field Bq. The Hamiltonian of that system 
in the rest frame is 



(13) 



where /i is the magnetic moment of the particle and 
is the Pauli matrix 



1 

-1 



Now let 



l*(0)) 



:os(6l/2)\ 
dn(0/2)) 



(14) 



(15) 



The generic derivation of the AA-phase for this system 
is given in numerous papers and books (see [Tsj ) and 
yields the famous result 7 — ±7r(l — cos(0)), the total 
phase (j) = ±7r being found from physical considerations. 



Alternatively since Aijt (and 
is known, one can calculate r 

Tft-i(«'(0)|iJ|*(0)) = -7rcos(6'). Moreover 



in fact A) 
from Eq.®, 
so that 



Eq. 



entails that the only possible total phase is solution 
oi (j) = 27r(n — 5) = 27r(m -(- i) with n and m two 
different integers, that is = ±7r. We thus recover 
7 7r(l — cos(0)). 

It is also possible to calculate the AA-phase from 
only one element of A$ using the fact that (f) = ±7r 
is known from physical considerations. We know that 
(1'(0)|ff|5'(0)) = -^Socos(6'). If we choose the eigen- 
value A+ = /iSo, Eg. pTj) gives us that 7 = 0(1-1- 
cos{6)) = ±7r(l -I- cos{9)) = ±7r =F Trcos{0) = q=7r =F 
TTcos{9) + 2n = ±7r(l — cos{d)) as it should. If instead 
we use, A_ = —fiBo, we have likewise immediately that 
7 ±7r(l - cos(6l)). 



B. Inflnite-size arbitrary Hamiltonian 

Consider now an infinite Hamiltonian and a state j^*) 
so that we know two different non-zero elements Ai and 
A2 of A^. We further assume that this Hamiltonian 
yields a cyclic motion of the state jvP). Our goal is to 
show how one can derive the geometric phase and con- 
strain the possible total phases, periods and the unknown 
part of the spectrum using these two eigenvalues. Fur- 
ther informations such as the knowledge of the period 
from physical considerations or experimental results al- 
low to refine these constraints. The Hamiltonian is given 
by the infinite matrix po| 



/Ai 

A2 

jj^ Bi B2 

Ci C2 

V : ; ; ; 



The geometric phase is given by Eq.© 



(16) 



7 = <^+^(i/) 



(17) 



where (H) indicates the expectation value of H calcu- 
lated for any state |^'(i)) located on the closed loop 
formed by the cyclic time-evolution in the projective 
Hilbert space. 

Consider first the most general case where none of the 
total phase and the period are known. Using Eq.®, we 
can constrain through the relations 



2irn — Ai — = 27rm — A2 
n 



(18) 
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Using the first relation, we find t as a function of (j) which 
allows to find (j) by introducing the expression of t in the 
second relation. We find 



27r 



T = 2T:h 



Aim — 
Ai - A2 ' 

n — m 

Ai ~ A2 ■ 



(19) 
(20) 



Now let us change of gauge to get = and thus 
Aim — K2n. This does not change the period nor the 
geometric phase and only shifts the whole spectrum of 
the Hamiltonian without affecting the dynamics of the 
system. The period is now expressed as 



T ^ 2Trh—. 
Ai 



This gives the geometric phase 



7 = 27T—{H)[2n], 
7 = 2^^(i/)[2^]. 



(21) 

(22) 
(23) 



As the geometric phase is defined modulo 2tt, if (if)/Ai.2 
is rational pTj . there is only a finite number of values of 
n and m that will yield different results up to 27r phase 
factors. It is interesting to remark that Eqs.® and (fT9|) 
imply that A1/A2 is a rational number, which, as we 
know from Sec lII B H and Sec lIIB 3l is consistent with the 
assumed fact that the system undergoes a cyclic motion. 

We now turn to the unknown part of the spectrum. 
Any eigenvalue of Aq, must fulfill and equation similar to 
Eg. psp and can thus be constrained as is now known 
to be 0. This leads to 



Aku = Aik' 



(24) 



k' being the integer entering Eq. (jl8p when written for A^. 
Such a procedure can be carried using other eigenvalues 
if some are known beyond Ai and A2. This will lead to 
new relations on the unknown A^ , (jj and r which further 
constrain their possible values. 

This illustrates that it is possible to partially know the 
geometric phase, (j), t and the unknown part of A^ from 



an incomplete spectrum from the condition of periodic- 
ity of the considered state \^p). If a total knowledge of 
the geometric phase is required, it is also possible to di- 
rectly impose constrains on the spectrum of H, as we 
now illustrate. 



C. Fabry-Perot Cavity with a Movable Mirror 

In this section we consider the problem of a single mode 
of the light field inside a Fabry-Perot cavity with one 
movable mirror. This system has been studied both the- 
oretically and experimentally in great detail in the con- 
text of the emerging field of cavity optomechanics. It is 
characterized by a large variety of quantum mechanical 
features, such as non classical states [ij] or entanglement 

m, [n, m. 

This system is described to an excellent approximation 
by the Hamiltonian j25| 

+ huj„,h% - hga)a{b + 6^), (25) 

where at (a) is the bosonic creation (annihilation) 
operator for the cavity field mode. Similarly, the 
mirror is treated as a quantum harmonic oscillator of 
frequency a;,„ with 6^ (&) the bosonic creation (anni- 
hilation) operator, and g is the opto-mechanical coupling. 

We first compute the geometric phase in this system 
via the generic method. The evolution operator U (t) was 
calculated and disentangled in Ref.[l3] as 

U{t) ~ cxp[— ira^ aujriit]cxp[ik'^ {a'^ a)'^ {u!,nt ~ sin u!,nt)] 
xcxp[ka'^a{rib^ — ri*b)]cxp[^ib''buijnt], 



where r = ujf/ ujm , = 5/ and 77 = 1 — e 
sider the initial state 



Con- 



(26) 



n=0 



where \n) f is a Fock state of the cavity field and |/3)m 
a coherent state of the mirror. Following the state 



|*(t)) at time t is 



00 



'2 

n=0 



r 



As in Ref.[l9[, we assume that r is an integer and we 
choose k so that fc^r = 2qTr and r = 2p7r, where p and 
q are the smallest integers satisfying fc^ = q/p (this is 
possible by constraining for example the length L of the 
cavity) . In that case the system's motion is periodic with 



total phase 4> = 27r. 



We still have to calculate d\'i'{t))/dt from Eq.([27I) and 
compute Eq.Q. The resulting geometric phase, calcu- 
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latcd in Ref. [T^ , is 



7 = 27r 
with (h) f 



^(r-2fci?e(/3))(n)/ + ^ 



(27) 



ErT^O^ICnl 



We now rederive that same result using our new 
method. We remark that G A* is an eigenvalue of H 
with eigenvector |0) field |0) minor = |00), thus (j) = 2it. Fur- 
thermore, using the initial state of Eq. ((26|) . the Hamilto- 
nian of Eq. ([25|) allows us to compute easily 



(f( 



|i7|*(0)) = nw„, [{r ~ 2kRe{(3)){n) f + 



and therefore the most general expression for the geo- 
metric phase is 



7 = 27r + u}mT{r — 2kRe{f3))(n) f + cj„ 



(28) 



The remaining difficulty is to find the period r. But we 
can remark also that any state of the form |0n), where 
n is a positive integer, is eigenstate of H with eigenvalue 
hwrnti G A*. Thus Eq.® together with (f) ^ 2tt gives us 
directly that fiwmLCM (^AE^^^ = p is a non-zero integer. 



Thus T 
Eq.dMl) 



27ih 



P 



2tt- 



Therefore we obtain from 



7 



27r 



^(r-2fci?e(/3))(n)/ + ^ 



(29) 



This expression holds without any assumption on the pa- 
rameters. However we want to check if we can make ap- 
proximations similar to those performed in the generic 
method, so that we get Eq. (P7|) and Eq. ip^]) equal. 



First we calculate fc^r = '^''^ 2mLV P- ^^P" 
pose that r is an integer. As in the generic method, 

is a non-zero integer. 



we constrain L such that — A a 
This way, we see that there exists a non-zero integer q 
such that = q/p and by definition of t, these inte- 
gers are the smallest satisfying this equation. Thus with 
the assumptions that r is an integer and that L is con- 
straint so that 7^ — r% i is a non-zero integer, we obtain 
that Eq.(l27| and Eq.dlSl) are equal. 



D. 



Three-Mirror Fabry-Perot Cavity with One 
Movable Mirror 



The last example that we consider is a three-mirror 
cavity with a slightly transmissive and movable middle 
mirror of mass m and natural vibration frequency ajm- 
This system has been discussed in Refs.[26|. [27j and [1^. 



1. Generic Method 

Ref. [131 identifies a regime of coupling between the 
middle mirror and the left and right cavity fields that 
is linear in the position of the mirror for one field and 
quadratic for the other. Expressing the momentum p 
and the position q of the mirror in terms of bosonic 

and c as p = 

f c) the Hamil- 



creation and annihilation operators 
i\/2nihw„i{c^ — c) and q = 



2mu. 



(ct 



tonian is 



H = hjjDO^a + hwsh^h + hCDa^a{c + c^) + {huj,m + hCsh'^h)c^c + ^^h^h{\ + + c^-^). 



(30) 



r 



where (a) and fet arc the bosonic creation (anni- 
hilation) operators of the light fields in, respectively, the 
left and the right cavity, and 



Cd = 

Cs = 



D 



^ 2mujm I h ' 

niLUm ' 



(31) 



with and ^5 constants, whose explicit forms are given 
in Ref.llal. We further introduce 



X = [u:^{uo^ + 2Csh%)\"\ 

S = {ujd — CD)a^o, + LUsb^b, 

ly = Coa^ax^^ [t^mX^^ (cosx^- 1) -isinxt] , 



\k\ = \smh-^{Csb^bx-'smxt)\ 



(32) 



Note that using realistic parameters for the system (see 
[Tsl). one obtain (x) — uJm- 

Using a semi-classical approximation for the light 
fields, the evolution operator corresponding to Eq. ipPl) 
was disentangled in Ref. [iBj . A full quantum treatment 
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leads to a similar evolution-operator 

U{t) = e'^^D{v)R{<^)S{K) (33) 

where D{v) = exp[z^c^ — v*c\ is the displacement oper- 
ator, i?($) = exp(i-j(ctc -f- cc^)) is a rotation operator, 
and S{k) = exp[K*c^/2 — kc^' V2] is a squeezing opera- 
tor. It has been shown in Ref. [l5j that the rotation due 
to R exactly oppose that due to S, so that Eq. ([33|l is 
equivalent to 

U{t) = e'^'D{v)S{\K\) (34) 

Now using the expressions of 5, v and k from Eg. ([5^ . 
one sees that \J(t) is periodic with period r ~ For 
that value of r, [/(r) = [/(O), so that </) 2ii^ This 
shows the difhculty of working with the generic method. 
Indeed, obtaining the geometric phase directly from [/(i) 
leads to rather involved expressions for both |^'(t)) and 

2. Obtaining the geometric phase 

The results for t and obtained in the previous section 
can also be obtained from H using a simple evaluation 



Eq. ([36|l is further simplified if we consider the initial state 
to be a product of coherent states |^'(0)) = \a)a\P)b\lJ)nn 



7 = 27r 



IV. CONCLUSION 

We have shown how, for time independent Hamiltoni- 
ans, it is possible to derive the geometric phase arising 
from a cyclic motion of the system without calculating 
the evolution operator. We have proposed several ways 
to compute t, and the resulting geometric phase 7 re- 
spectively through Eq.®, Eq.®, Eq.® and Eq. pT]) . 
This led to the conditions for an Hamiltonian to yield a 



of some of its eigenvalues, i.e. without calculating U{t). 
Consider for example an initial product state 



|*(0)) = ^"1")'^ Mn\n),n- (35) 

n— n— n— 

First, we remark that the state 

|0)fioid a|0)fioid b|0)mirror = |000) is an dgcnstate of 
H with eigenvalue G A^. Thus we have </) = 27r. 
Furthermore the state |00n) (where n is a positive 
integer) is eigenstate of H with eigenvalue hujmn € A*. 
Therefore there exists a non-zero integer p such that 
T = 2tt-^. This gives the period of cyclic motion with- 
out any assumption. It is nevertheless possible to find 
specific values of p using approximations. For example, 
considering that {2Csb^b) ^ ujm and {Cdo^o) ^ 'j-'m, 
(jjo/uim and LOs/ijJm integers, yields p = 1 (free field of 
frequency uj„i). 

To calculate the geometric phase, it remains to intro- 
duce T and (f) in Eq.® and to choose an initial state. The 
expectation value of at t = is easily calculated from 
Eq. ((30l) and Eq. pS]) , and the most general expression for 
the geometric phase given by Eq.® is 



(36) 



in which case 



(37) 

I 

cyclic motion of a given state. We have further shown 
how using the periodicity condition and knowing only a 
part of the spectrum of the Hamiltonian it is possible to 
derive the geometric phase, with the remaining uncer- 
tainty contained in the total phase or the period of the 
cyclic motion. This also permits to build constraints on 
the unknown part of the spectrum of H. Through the 
spin — ^ particle and the two mirror vibrating Fabry- 
Perot cavity examples, we have shown that our method 



= 27r 



l+p['^{d) + ^{b)+ 2^(5) V V^i?e(Af„M„+i) 

\ LOm OJm OJm — ^ 



(m) + — (fe)(m 



Cs 



ri=0 



(fo) 1 + 2^ V(n + 2)(n + l)i?e(M„A^„+2) 



n=0 



I 

\+p\a\^ f ^ + 2^i?e(M)) +pW(^ + + 2i?e(M)^)) 

\UJrn / ^ / 
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gives the same result as the generic method. 
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APPENDIX A: EXPRESSIONS OF r AND 

In this Appendix we derive the expressions for r and 
(j) respectively given in Eq.([6]) and Eq.®. 

Let B — {\4>k)} be a basis in which H is diagonal 
and A = {A^} the corresponding set of its eigenvalues. 
Let C B the be smallest set of eigenvectors needed 
to decompose a state \'^) on B and let C A be the 
corresponding set of eigenvalues. Finally, let AE^ be 
the set of non-zero energy spacings in A^, i.e. AE>j, = 

In the basis B, U{t) is diagonal and the elements of 
the diagonal have the form exp[— i/i~^At]. By definition 
of the period, r is given by 

T = mintjfo ^ 0, C/(io)|*) = e^^C/(0)|«')}, (Al) 

which, in the diagonal basis, means that there exists 
such as for all the elements A G Aq, , there exists an integer 
n such that 

h'^Xta = -(j) + 2m: (A2) 

and T = mintgjto}- Let us say that the above condition 
has at least one non-zero solution, i.e. that the system 



effectively undergoes a cyclic evolution with total phase 
<j>. We now subtract Eq. (|A2p for two different elements of 
A* (we suppose that A^ has at least two different eigen- 
values, otherwise the state does not evolve with time). 
This leads to the following condition that is true for any 
pairwise different eigenvalues A^ and Ai elements of A$ 

3m eZ such that h^^toAEkj = 2m7r, (A3) 

one can see that any to/(27rft) solution of Eq. (jA2p must 
be proportional to a multiple of any AE^^ and there- 
fore be among the common multiples of all elements of 
AE^^ . T I {2T:h) being the smallest quantity to fulfill this 
condition, t is proportional to the least common multiple 
(LCM) of AE^\ We obtain Eq.® given in Sec HTBl] 

T = 2TTh LCM (AE^^) (A4) 

Note that in the case where is a integer multiple 
of 2tt, Eq.jAll reduces r = 2nh LCM(A-i), A e A*, 
as Ea. (|A2[) shows that for all A e A;^, there exists an 
integer m such that to = 2rm:hX^^ . Similarly, in the 
case where all the elements of A^ are equal (say to 
A), T simply reduces to 2'Kh\^^ (period of the complex 
exponential ex^[—ihr^Xt\. 

As T is known, we can derive the expression of the total 
phase 4) inserting Eq.(|^4| in Ea. (|M)) : 

(/. = 27r[n-ALCM(A£;-i)] (A5) 

where n is an integer and depend on the considered eigen- 
value A G Aq,. Note that this expression is valid for any 
eigenvalue element of A.p . 
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